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Abstract 



In the classic rumor spreading model, also known as susceptible-infected- 

■ removed (SIR) model in which rumor spread out of each node proportionally 
QQ ■ to its degree and all edges have uniform tie strength between nodes. The ru- 

i mor spread phenomenon is similar as epidemic spread, in which all the informed 

\ nodes spread rumor by informing their neighbor nodes. In earlier rumor spreading 

models, at each time step nodes contact all of their neighbors. In more realis- 
tic scenario it is possible that a node may contact only some of its neighbors 
^ I to spread the rumor. Therefore it is must in real world complex networks, the 

■ classic rumor spreading model need to be modified to consider the dependence 
of rumor spread rate on the degree of the spreader and the informed nodes. We 
have given a modified rumor spreading model to accommodate these facts. This 
new model, has been studied for rumor spreading in complex networks in this 
work. Rumor spreads nonlinearly by pairwise contacts between nodes in complex 
networks at degree dependent spreading rate. Nonlinear rumor spread exponent 
a and degree dependent tie strength of nodes affect the rumor threshold. The 
de gree exponent tie strength between two nodes is (/cj/cj)^, where ki and kj are 
degrees of node i and j and /? is tie strength exponent. By using the above two 
exponents in any complex network gives rumor threshold as some finite value. In 
the present work, the modified rumor spreading model has been studied in scale 
free networks. It has been found to have a greater threshold than with previous 
classic rumor spreading model. It is also found that if a and (3 parameters are 
tuned to appropriate value, the rumor threshold becomes independent of network 
size. In any social network, rumors can spread may have undesirable effect. One 
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of the possible solutions to control rumor spread, is to inoculate a certain frac- 
tion of nodes against rumors. The inoculation can be done randomly or in a 
targeted fashion. We have used modified rumor spreading model over scale free 
networks to investigate the efficacy of inoculation. Random and targeted inoc- 
ulation schemes have been applied. It has been observed that rumor threshold 
in random inoculation scheme is greater than the rumor threshold in the model 
without any inoculation scheme. But random inoculation is not that much effec- 
tive. The rumor threshold in targeted inoculation is very high than the rumor 
threshold in the random inoculation in suppressing the rumor. It has also been 
observed that targeted inoculation scheme is more successful to stop the rumor in 
this model based on scale free networks. The proposed hypothesis is also verified 
by simulation results. 

Keywords : Complex networks; Scale free networks; rumor spreading; nonlinear 
rumor spread; tie strength in complex networks; random and targeted inoculations 



1 Introduction 



Many researchers have discussed, how the properties of networks affect the dynamical 
process taking place in networks. In the recent years, complex network structures 
and their dynamics have been studied intensively [U IH [THl [TBI [191 EOl EEl [29]. By 
analyzing different real world networks e.g. Internet, the www, social network and so 
on, researchers have identified different topological characteristics of complex network 
such as the small world phenomenon and scale free property. An interesting dynamical 
process in complex networks is the epidemic spreading. Work on epidemic spreading has 
been done by many researchers [171 [I9]. There are two models to define the epidemic 
spreading. In first model- two state susceptible-infected-susceptible (SIS) model, a 
susceptible node can become infected and an infected node can recover and return 
to susceptible state, e.g. computer viruses [6l [2l]. The second model, is susceptible- 
infected-recovered (SIR) model, it defines that the infected nodes will become dead 
pn [T9] . It is different from SIS model in the fact that the infected nodes will not return 
to the susceptible state but can attain recovered status. The epidemic spreading on 
complex networks has been discussed with the other models also [23]. Rumor spreading 
is also similar in nature as epidemic spread. 

Rumors are information that circulate without officially publicized confirmation i.e. 
without has official reputation. Rumors are being used as a special weapon of public 
opinion and can pose a tremendous impact on social life. Certain rumors can affect 
the social stability seriously. In order to improve the resistance of the community 
against undesirable rumors, it is essential to develop deep understanding of the mech- 
anism, underlying laws involved in rumor spreading, establishment of an appropriate 
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prevention and control system to generate social stability. Sudbury studied, first time 
on the spread of rumors based on SIR model [28]. Another standard model of rumor 
spreading, was introduced many years ago by Daley and Kendal [9]. Its variant was 
introduced by Maki-Thomsan [Mj. In Daley-Kendal (DK) model homogeneous pop- 
ulation is subdivided into three groups: ignorant (who don't know rumor), spreaders 
(who know rumor) and stifier (know rumor but do not want to spread it). The rumor is 
propagated through the population by pairwise contacts between spreaders and other 
individuals in the population. Any spreader involved in a pairwise meeting attempts 
to infect other individual with the rumor. In case this other individual is an ignorant, 
it becomes a spreader. If other individual is a spreader or stifier, it finds that rumor is 
known and decide not to spread rumor anymore, thereby turning into stifier. In Maki 
Thomsan (MK) model when spreader contacts another spreader, only the initiating 
spreader becomes a stifier. DK and MK models have an important shortcoming that 
they do not take into account the topology of the underlying social interconnection 
networks along which rumors spread. These models are restricted in explaining real 
world scenario for rumor spreading. By considering the topology of network, rumor 
model on small world network |18l [31], [32] and scale free networks have been de- 
fined. Therefore, as long as one knows the structure of spreading networks, he can 
figure out variables and observable to conduct quantitative analysis, forecasting and 
control of rumor spreading. The most important conclusion of classical propagation 
theory is existence of critical point of rumor transmission intensity. When an actual 
intensity is greater than critical value, the rumors can spread in networks and persis- 
tently exist. When the actual intensity is less than the critical value, rumors decay at 
an exponential rate and this critical value is called rumor threshold. In previous studies 
of rumor spreading, underlying network topologies were scale free which is also observed 
in many real networks. But some more properties of real world networks have not been 
considered. In the earlier rumor spreading models, rumor transmission rate was fixed. 
In the real scenario, it will be different among different nodes. It will depend on tie 
strength between nodes and node's strength. In Internet, weight implies the knowledge 
of its traffic flow or the bandwidth of routers [22], in the world wide airport networks 
it can deflne the importance of an airport [S| and so on. In case of rumor spreading, 
the strength can indicate the frequency of the contact between two nodes in scale free 
networks. Greater the strength, the more intensely the two nodes are communicating. 
Chances of spreading rumors tend to differ among individuals and laws of spreading 
in social networks with different topologies are also different. Each informed node can 
make contacts with all of its neighbors in a single time step. In other words we can say 
that each informed to number of nodes, can spread information to nodes equal to its 
degree. In real case, an informed node can't make contact to all of its neighbors in single 
time step. Studies on small world networks found that compared with regular network, 
small world network has smaller transmission threshold and faster dissemination. Even 
at small spreading rates, rumors can exist for long. Studies on infinite-size scale free 
networks have also revealed that no matter how small transmission intensity be, the 
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rumors can be persistent as positive critical threshold does not exist. [25| 126] . 

In previous studies on rumor spreading in scale free networks, it has been assumed 
that larger the nodal degree, the greater the rumor spreading from the informed node, 
i.e. the rumor spread is proportional to the nodal degree. With these assumptions 
for SIR model, for scale free networks with sufficiently large size, the rumor threshold 
Ac = can be zero. Yan et al. have demonstrated that the asymmetry of infection 
plays an important role [30]. They redistribute the asymmetry to balance the degree 
heterogeneity of the network and found finite value of epidemic threshold. Zhou et al. 
concluded that this hypothesis is not always correct [2[ |33] . In rumor spreading contact 
networks, the hub nodes have many acquaintances; however they cannot contact all 
their acquaintances in single time step. They assumed that the rumor spreadness is not 
equal to the degree but identical for all nodes of the scale free networks and obtained the 
threshold Xc = \, where A is the constant infectivity of each node and is not equal to the 
degree of node. Recently, Fu et al. [10] have defined piecewise linear infectivity. They 
suggested if the degree k, of a node is small, its infectivity is ak otherwise its infectivity 
is a saturated value A when k is beyond a constant A/a. In both constant and piecewise 
linear infectivity, the heterogeneous infectivity of the nodes due to different degrees has 
not been considered. While in scale free networks heterogenity in nodal degree is very 
common. There may be nodes with different degrees, which have the same infectivity, 
and there will be a large number of such nodes if infectivity does not saturate or the 
size of network is infinite. 

In scale-free networks, a small number of nodes have very high degrees. Although, 
random immunization strategy works very well in homogeneous random networks, but 
this strategy is not effective in preventing a rumor in scale free networks. Hence a 
new immunization strategy needs to be developed which is able to recover the rumor 
threshold. One of the most efficient approach is to immune the highest degrees nodes, 
or, more specifically, to immune those nodes (hereafter termed as hubs or hub nodes) 
which have degrees higher than a preset cut-off value k^ . Such a strategy is known as 
targeted immunization [3 d El [IH [251 IMl EZ] . 

In order to control the epidemic spread of rumors, inoculating the nodes is an option. 
Random inoculation was found to be ineffective for all the complex networks (includes 
scale free networks) by Pastor-Satorras and Vespignani [26]. On the other hand, tar- 
geted inoculation is successful in arresting the rumor spread in scale free networks [26j . 
Acquaintance inoculation method proposed by Cohen et.al. [13] was found effective 
for both scale free networks as well as networks with bimodal degree distribution. Site 
percolation is another process inversely similar to inoculation [19]. Inoculated node can 
be considered to equivalent of removed site in site percolation. After the percolation 
threshold, the network will have maximum size of giant component [IH]. Thus if the 
network operates below the percolation threshold, the rumor spread will not happen as 
outbreak. In other worlds, inoculation ensures the operation of network below perco- 
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lation threshold. Random inoculation usually requires inoculation of large number of 
nodes for being effective. If nodes with higher connectivity are targeted for inoculation, 
the same effectiveness can be achieved with smaller number of inoculated nodes. But 
it needs knowledge of nodes which have higher connectivity \13\ . 

In this work, we have investigated rumor spread for the scale free network while 
considering varying tie strengths between nodes. Further we have assumed that a non 
linearly varying number of neighbors are infected with the rumor in each time step. 
While in the earlier models |T71 118], tie strength has been considered to be uniform, 
and a constant number of neighbors have been assumed to be infected in each time step 
by each node. If a node has K neighbors, in the earlier models, in each time step, all the 
K neighbors will be infected. We have modified the earlier SIR model given by Nevokee 
|18l| and included a rumor spreading exponent a . In this work, K"' neighboring nodes 
will be infected in each time step. Here a is the spreading exponent where, < a < 1. 
We have used Barabasi- Albert (BA) model [1] to create scale free networks with power 
law distribution of nodal degree, and then used the proposed strategy in them to study 
the rumor spread. Scale free networks have been specifically chosen as they are much 
more heterogeneous than the small world or the random network models, thus forming 
a good candidate for testing our proposition. 

The dynamical differential equations have been used to represent the modified model 
for information spread. The equations have been used to find the threshold and study 
the rumor propagation behavior. The result have also been verified by the simulations 
By choosing the appropriate values of a and /3, finite non zero threshold value can be 
found for the scale free networks. It is found that the a is more sensitive than the /3 for 
spreading rumor in a large scale free network. Finally, the rumor threshold has been 
calculated after applying the random and targeted inoculation of nodes to suppress the 
harmful rumor in the scale free networks. In the targeted inoculation scheme, the rumor 
threshold has been found to be larger. The targeted inoculation has also been effective 
in the scale free networks to suppress rumors. 



2 Classical SIR model 



Classical SIR model is one of the most investigated rumor spreading models for complex 
networks. In this model, nodes are in one of the three categories — ignorants (the nodes 
who are ignorant of the rumor), spreaders (those who hear the rumor and also actively 
spread it) and stifiers (the nodes who hear the rumor but do not spread it further). The 
rumor is propagated through the nodes by pairwise contacts between the spreaders and 
other nodes in the network. Following the law of mass action, the spreading process 
evolves with direct contact of the spreaders with others in the population. These con- 
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tacts can only take place along the edges of undirected graph of complex network. If 
the other node is the spreader or stifler then the initiating spreader becomes the stifler. 
The classical SIR model has been studied by M. Nevokee et al. [161 IlH] for heterogeneous 
population (nodes have different degrees). In this paper I{k,t), S{k,t), R{k,t) are the 
expected values of ignorants, spreaders and stifler nodes in network with degree k at 
time t. Above rumor spreading process can be summarized by following set of pairwise 
interactions. 



h + S2^Si + S2, 

(when spreader meets with the ignorant, it becomes spreader at rate A) 
Si + R2 ^ Ri + R2, 

(when a spreader contacts with stifler, the spreader 
becomes a stifler at the rate 6) 

+ 5*2 ^ -Ri + 5*2, 

(when a spreader contacts with another spreader, initiating spreader 

becomes a stifler at the rate 6) 

S^R. 

{a is the rate to stop spreading of a rumor spontaneously.) 

Let p\k,t) = I{k,t)/N{k),p'{k,t) = S{k,t)/N{k), p''{k,t) = R{k,t)/N{k) are the frac- 
tion with respect to total nodes, of ignorant, spreaders and stifler nodes respectively 
with degree k at time t. These fractions of the nodes satisfy the normalization condition 
p\k,t) + p'^{k,t) + p^{k,t) = 1 where, N{k) represents the total number of nodes with 
degree k, in the network. Nevokee et al. [18] proposed the formulation of this model 
for analyzing complex networks as interacting Markov chains. They used the frame- 
work to derive from the first-principles, mean-field equations for the dynamics of rumor 
spreading in complex networks with arbitrary correlations. These are given below. 

= -kXp^ik, t) J2 Pmp'{l. t); (1) 

dp'{k,t) 

Jt 

/ I 
+p^{l,t))Pm-^p'ik,t); (2) 
dp''{k,t) 



k\p\k, t) pmp^i, t) - ksp^k, t) J2 (p'(^ t) 



dt 



k5p^{k, t) 5^(p^(/, t) + p\l, t))Pm + ap\k, t)- (3) 



Where conditional probability P{l/k) is the degree-degree correlation function that 
a randomly chosen edge emanating from a node of degree k leads to a node of degree 
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/. Here, it has been assumed that the degree of nodes in the whole network are uncor- 
related. Therefore, degree-degree correlation is P{l\k) = ^-j^ where P{1) is the degree 
distribution and (k) is the average degree of the network. Nevokee et al. have shown 
that the critical threshold for rumor spreading is independent of the stifling mechanism. 
The critical threshold found by him was Ac = It is same as found for SIR model 
[El [17]. Hence, it implies epidemic threshold is absent in large size scale free networks 
{{k"^) — J- oo, Ac — )■ 0). This result is not good for epidemic control, since the epidemics 
will exist in the real networks for any non zero value of spreading rate A. 

The "scale-free network" (Internet, www etc.) [Ull] describes the degree distribution 
as a power law, i.e., the probability of a node having degree k 

P{k) oc k-"' 

Where, 2 < 7 < 3. Degree distribution is plotted on a log - log scale where, 7 is the 
slope of the degree distribution, which is a straight hne. One can construct a scale- 
free network by starting with a few nodes moand then by adding a new node with 
m {m < mo) edges per node at each time step. The probability H^, of node v being 
connected to a new node depends on its degree j„. The higher its degree, the greater the 
probability (H„ = j^/ Ylikik) that it will get new connections (preferential attachment). 



3 Modified rumor spreading model 



The real world networks can have the intimacy, confidence etc. between the nodes. Un- 
like previous studies where each node can spread the rumor with constant transmission 
rate A. In this study, we have concentrated the rumor spreading model with nonlinear 
rumor spread and transmission rate between two connected nodes is a function of their 
degrees. Here, to remove the complexity stifling rate 5 is being considered 0. Based on 
this assumption, we define, 

'^-^^ = ~kp\Kt)Y,P{l\k)p\l,t)^\,,- (4) 
'^^^^'''^^ =-ap^ik,t) + k{l-p^ik,t)~p^ik,t))J2Pmp'il,t)^\i,; (5) 



dt 

dp''{k,t) 
Jt 



ap%k,t). (6) 



where Q{k) represents the rumor spreadness of a node with degree / and Xik represents 
the rumor spreading rate from nodes of degree / to nodes with degree k. 
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3.1 Tie strength in complex networks 



The topological properties of a graph are fully encoded in its adjacency matrix A , whose 
elements aij {i 7^ j) are 1 if a link connects node i to node j, and otherwise. The 
indices i, j run from 1 to N, where N is the size of the network. Similarly, a weighted 
network is entirely described by a matrix W whose entry Wij gives the weight on the 
edge connecting the vertices i and j {wij = 0, if the nodes i and j are not connected). 
In this study we will consider only the case of symmetric weights {wij = wji) while the 
undirected case of the network is considered [5]. In a call network if two nodes call 
each other for a long duration then weight of the connecting edge will be high and it 
shows high tie strength between them [21]. Here weight of the edge in terms of total 
call duration defines the tie strength between the nodes. It has also been observed 
in the dependence of the edge weight Wij to define strength between nodes with end 
point degrees ki and kj. Weight as a function of the end-point degrees can be well 
approximated by a power-law dependence, 

Wij = h{kikjY 

where, /3 is the degree influenced real exponent which depends on the type of complex 
networks and 6 is a positive quantity. When /3 > then rumor transmit to high degree 
nodes and when /3 < then rumor will transmit to low degree nodes. Further, if /3 = 
there will be degree independent transmission. 

It has been observed that the individual edge weight doesn't provide clear view of 
network's complexity. A detailed measurement if tie strength using the actual weights 
is obtained by enhancing the property of a vertex degree fcj = "Y^- a^j in terms of 

the vertex strength Si = X^jLi '^jj'^u (total weights of their neighbors). Therefore, 
there is a coupling between interaction strengths of the nodes with the counterintuitive 
consequence that social networks are robust enough to the removal of the strong ties 
but fall apart after a phase transition if the weak ties are removed [^5.j. Therefore, we 
can measure the strength of a node of degree k for scale free network, 

Sk = k^P{l\k)wki 
IP{1). 



{k) 



-Wkl 



*^(*'-) (7) 



Here, it has been considered that the rumor spreading model, where rumor transmis- 
sion rate in contact process between a spreader node and an ignorant node influenced 
by their degrees. If Wki is the tie strength between ^-degree node and /-degree node for 
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{k,l) edge, Sk is the node strength with degree k. In scale free network each node of 
degree k, there is a constant rumor transmission rate Xk. Therefore, rumor transmission 
rate from k-degree node to /-degree node is given by the proportion of Wki to Sk- Hence, 
Xki can be defined as, 

Xki = Afc^ (8) 

We can see in Eq. ([S]) that by increasing the proportion of Wki/Sk, the more possibihty 
of rumor transmission rate can be increases through the edge. In the present work, 
uncorrelated networks have been considered, hence A^; = Xl^ (k) / (k^^^). In this model, 
rumor spreadness, Q{k) = fc" where < a < 1, it defines that each spreader node may 
contact with fc" neighbors within one time step. Therefore, spreadness of a rumor will 
vary nonlinearly with the growing degree k. In Eqs. by putting a = 1 (without 

the lose of generality), we can write rumor equation for Q{k) and Xik can be written as, 

^^ = -^p^(^>^)E^"^(Op^(M); (9) 

= -P'ik, t) + j^^P'ik, t) ^ rP(/)p^(/, t); (10) 

dp''{k,t) 



dt 



p\k,t) (11) 



Here, ^(t) = Xlfc ^) has been used as auxiliary function. After solving 

rumor Eqs. (fTOl)-(fTT!) with initial conditions p*(A;,0) ~ l,p^(A;,0) ~ 0,p^(A;,0) ~ we 
get. 



p\k,t)=e^^ (12) 
where, 

^(t)= / $(t)dt = ^A;°P(A;)p"(A;,t) (13) 
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4 Rumor threshold of the modified model 



Time derivative of Eq. f|T3|) can be obtained by, 
dt ^ ^ ' dt 

k 

from Eq. ([12]), 

k 
k 

= J2 ^"^(^) - E k"Pik)p'{k, t)-J2 k^Pik)p\k, t)) 

k k k 

= {k"") -"^{t) -J2k"P{k)p\k,t)) 

k 

= (fc") -^(t) -^A;°P(/t)e ^ ^ (14) 

k 

In the infinite time limit, i.e., at the end of rumor spreading, we will have S'fc(oo)=0, 
limt^oo^{t) — )■ Tp and limt^aod'^ /dt=0, 

= (A;") - ^ - ^ fc"P(A;)e 

k 

^ = (fc") - ^ fc"P(fc)e (15) 

A; 

We can conclude that \& = will always be the solution of above equation. It 
has been found that right hand side of Eq. f|T5|) is a convex function. Therefore the 
derivative of Eq. ( |T5i) with respect to \E' will always be in order to achieve greater than 
one to get non zero solution, 

(('^°) - Yl k''P{k)e <^ )|*=o > 1 (16) 



Hence, 



i:mA(^ = Al^>l (17) 



Using above equation, the rumor threshold can be defined as. 



A. 
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It is interesting to note that, by putting, a = 1 and /3 = in Eq. (ITSl) . the threshold 
for this model reduces to {k)/{k'^) for classical rumor spreading mo del [26]. 



When, t — )■ oo spreader nodes will be 0, {p'^ik, oo) = 0) and from Eq. (fT2l) . p\k, oo) = 
g"(fei+^ . Therefore, final size of rumor i? at t — t- oo {limt^aoP^ {k^t) = R), 

R = Y.P{k)p\k,^) (19) 

k 

= 5^P(fc)(l-p^(A:,oo)) 

k 

= ^P(A;)(l-e^^^ ) 

= l-^P(A:)e^^^ (20) 



(21) 

In the real world complex networks rumor spreads on a finite size complex networks. 
It may be possible that size of scale free network is very large. The maximum or minimum 
degree of scale free network is mentioned by fcmax or ^mm- Pastor et al. found that the 
epidemic threshold Ac for k^ax for SIS model on bounded SF networks with P{k) ~ 
^-2-7 ^ < 7' < 1. They assumed that with the soft and hard cut-off kmin and kmax, 
when a = 1. The hard cut-off denotes that, a network does not possess any node with 
degree k > kmax- As kmax of a node is network age, defined in the terms of number of 
nodes N, 

kmax = kminN ^'+^ (22) 

The normalized degree distribution is defined by, 

-L yf^max I rimin) 

Where d{x) is a heavy side step function 



In modified rumor spreading model if a = 1 and /3 = then it will be classic rumor 
spreading model and degree distribution in scale free networks P{k) = k"^ where, 
2 < 7 < 3, therefore 

\c{k^ax) = (24) 

= (25) 

— 7 r)\] i,kinax / kmin)'^ (26) 
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Eq.( l22|) is modified for tlie given scale free network as, 

XJ{N) ^ , (iV)(^-3)/(7-i) (28) 

for 7 = 3, 

7c'(iV) ^2[A;™„/n(iV)]~i (29) 
Eqs. ([28])- ([29]) shows tliat A', -> if ^ oo 

In modified rumor spreading model, nonlinear rumor spread is considered with 

e{k) = r , 



(_^) g + /3 - 7 + 2 [{k^ajkminy-^^^ - 1] 

/3-7 + 2 [(fc^„.A™0"+^-^+' - 1] ^ ' 



Theorem 4.1 /n classic rumor spread model (a = 1, P = 0) threshold is smaller than 
the modified rumor spread model (0 < a < 1 and P 0). 

The proof is given in the next section after lemmas. 

Lemma 4.2 When the size of network (N) increases, the value of critical threshold 
Xf>0 for a + /3 + 2 < 7, otherwise it will approaches to 0. 

Proof Since kmax/kmin = N'^^"'~^, therefore kmax/kmin increases when A'" increases, it 
becomes infinity when -> oo. When, a + (3 + 2 < 7, (|nM-)/3-7+2 = ^kjna.\a+f}-^+2 ^ q_ 

The value of Xf will be positive. Here, /3 < (rumor transmission influenced to low 
degree nodes) is considered. Now from Eq. (130]) . we can conclude that Xf will be 
positive . For a + /3 + 2>7. A*— )-0 when N increases. It can be summarized as, 

{^min 7-^-2^ (kmax/kmin)'^ a/3 2^ q,_|_^_|_2>7 
4m^^^^5!^, a + /3 + 2<7 (31) 

kl^Z^ —r-r, — ^l T , a + /3 + 2 = 7 

Lemma 4.3 In given rumor spreading model when a + (3 + 2 < 7 then rumor spreading 
threshold A* is independent from the size of scale free network (N). 
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Proof It may also be defined using Eqs. (127I) - (13T!) in the term of the number of nodes 
N, 



A/(iV) = 4".?^^^^, a + /3 + 2<7 (32) 

K, '^min aln{N) ' u. ^ fj J 

Here, it is found that for a + (3 + 2 < •y, Xf is independent of N. 



Proof Now using lemmas (1) and (2)the theorem can be proved for a + P — j + 2 > j. 
The ratio of rumor threshold in classic model and given model is given as, 

Ac(iV) _ (2-7)(7-/3-2) ^33^ 



X*{N) (7 - 2)fcL-°)(a + /3 - 7 + 2)N(1 - y - (3 + 2)/(7 - 3) 

It has been found from Eq. fl33|) that ^"S^l, < 1 finite scale free networks. 

Therefore, it has been justified that rumor threshold Xf{N) is greater than the Xc{N) 
in finite size scale free networks. In finite size scale free networks, when < a < 1, 
/3 7^ and a + /3 + 2 > 7 then it is hard to spread rumor in comparison to networks 
have a = 1 and f3 = 0. Finite rumor threshold is possible for any size of networks as 
seen in Eq. ( 132|) . However, it will be when N approaches infinity. 



5 Random inoculation 



In random inoculation strategy, randomly selected node will be inoculated. This ap- 
proach inoculates a fraction of the nodes randomly, without any information of the 
network. Here variable g {0 < g < 1) defines the fraction of inoculative nodes. In the 
presence of random inoculation rumor spreading rate A reduced by a factor (1 — g). 
In mean field level, for the scale free networks in the case of random inoculation, the 
rumor equations are modified using initial conditions as, 

MM . -p.(M) + M!:!,.(M)tW. (35, 

dp''{k,t) 



dt 



p%k,t). (36) 
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Now \1/ is modified as (from Eq. (fT5|) ). 



^ = (fc") - ^ A;°P(A;)e <^^+'^> (37) 



Tfierefore, final size of informed nodes (R) is, 



-A(l-g)fcl+/^ ,., 

R=l-^P{k){l-g)e 0^'+") -g {31 

k 



To acfiieve nontrivial solution of \f' from Eq. (13 



A((fc«) _ ^ fc"P(A;)e"t^+r'*)|^^o > 1 (39) 

k 

Therefore, rumor spreading tfireshold in tlie case of random inoculation is obtained as, 

^4U) 

The relation between rumor spreading threshold, with inoculation (Ac) and without 
inoculation(Ac) can be defined as. 

It is to note that by applying random inoculation, the rumor spreading threshold (Ac) 
can be increased as seeen in Eq. (14T1) (i.e., Ac > Ac). 



6 Targeted inoculation 

Scale free networks permit efficient strategies and depend upon the hierarchy of nodes. 
It has been shown that SF networks shows robustness against random inoculation. It 
shows that the high fraction of inoculation of nodes can be resisted without loosing 
its global connectivity. But on the other hand SF networks are strongly affected by 
targeted inoculation of nodes. The SF network suffers an interesting reduction of its 
robustness to carry information. In targeted inoculation, the high degree nodes have 
been inoculated progressively, i.e more likely to spread the information. In SF networks, 
the robustness of the network decreases at the affect of a tiny fraction of inoculated 
individuals. 

Let us assume that fraction gk of nodes with degree k are successfully inoculated. 
An upper threshold of degree kt , such that all nodes with degree k > kt get inoculated. 
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Fraction gk of nodes with the degree k are successfully inoculated. The fraction of 
inoculated nodes given by, 




9k ^ if: (42) 



where < / < 1, and X]fe5'fe-P(^) — 9-i where g is the average inoculation fraction. 
Therefore, now rumor spreading equation is defined for targeted inoculation as, 

Ml -9.)»='t!,<(,, (43) 



dt 

Further, ^' for targeted inoculation. 



dt {k 

. . ^(1^,.(M)*W; (44) 

(45) 



A(l-9fc)fc^ +^ 



^(t) = (A;") -^A;"P(A;')e <^^+^> (46) 

k 

-J2k-P{k)e"'y^^)\^=o > 1 (47) 

k 

Therefore, rumor spreading threshold in the case of targeted inoculation is obtained as. 



Here, {gkk''+^+^)=g{k'^+^+^) + rjf, where rjf = {{gk - g)[{k''+^+^ - is the co- 

variance of gk and k'^'^^'^^. The cut-off degree k^ is large enough where rj/ < 0, but for 
small kc, gk — g and — have the same signs except for k's where gk — g 

and/or - is 0. 

Therefore rjf > for appropriate k^ 

Ac > (49) 

If fractions of inoculation in random and targeted are same then g — g, 

A~c > Ae (50) 

The above relation shows that in scale free networks targeted inoculation is more effec- 
tive than the random inoculation. 
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7 Numerical simulations: results and discussion 



The studies of uncorrelated networks have been performed using the degree distribution 
of scale free network. The size of the network is considered to be = 10^ and the degree 
exponent (7)=2.4. At the starting of rumor spreading, the spreaders are randomly 
chosen. In Fig. [H the final size of rumor R is plotted against rumor transmission rate 
for A^=100000, 1000 and 100 by tuning a and /3 as, 

• a + P = 0: Finite rumor threshold has been found. It has been observed that for 
different size of networks, constant threshold is there (after fixing the value of a 
and P). For the case a + /3 + 2 < 7, since 7 = 2.4. Therefore here it is interesting 
to see that finite threshold has been found which is independent from the size of 
network same as obtained from Eqs.f l3ip and (132|1 . 

• a + P = —1 : The simulation results are found same as above since a + /3 + 2 < 7 
with finite threshold and constant for any network size (for fix values of a and 

• a + P = 1: The results show some rumor threshold but approaches to as network 
size increases. For this case a + /3 + 2 > 7, since 7 = 2.4. Therefore, the threshold 
approaches to as network size increases. Similar results have been obtained by 
Eqs.(l3l])-(I32]). 

• a + P = 2: The simulation results are found same as above since a + /3 + 2 > 7 
with threshold approaches to as size of network increases. 

Using simulation results, of final size of rumor (i?), is plotted against time [t) in Fig. 
|2j It has been observed that rumor size increases exponentially as time increases and 
after some time it approaches a steady state, that will remain constant, since spreader 
density is at that time. It has also been observed that when a + /3 is low, then rumor 
size increases slowly initially but when a + P increases rumor size increases rapidly 
against time. While tuning a + P = —1 to a + P = 2 rumor increments faster initially 
(Fig. |2]). When ratio of a and P is high than the rumor size is also high. This result 
justifies that the a affects more final rumor size R than pit is seen from Eq. (fT8|) that 
when a is very small (0.1. 0.3) then the rumor threshold will be high. Than it is seen 
the final size of rumor will be too small when rumor transmission rate (A) is less than 
the rumor threshold (Ac). 

Critical rumor threshold is plotted against a in Fig. |3] while considering P = 0. 
Here, Ac decreased exponentially with the increase of a. It is maximum for a = 0.1 and 
almost at a = 1 for = 100000. Interestingly, this also happens in real life situation, 
when an informed node pass information to its maximum number of neighbors then 
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(a) R vs A for a + /3 + 2 > 7 



(b) R vs /3 for a + /3 + 2 < 7 



Figure 1: R vs X with a + /3 + 2 > 7 and a + /3 + 2 < 7 for different size of scale free 
networks 

rumor spreading will get outbreak in the network. However, the outbreak is hard to 
achieve when it passes rumor to less number of neighbors (~ 10-30%). Similarly, Ac 
has been studied against /9 at a = 1 in Fig. |3l It is found that /3 affects less rumor 
threshold for entire range except when /3 > 0. Further it approaches to 0. When size 
of the network (A^) increases then rumor threshold is decreased as shown in Fig. |3][a) 
and Fig. Mh)- 

In Fig. m final size of rumor has been plotted against /3 for a = 1, 7 = 2.4, 
3, and A = 1. It is seen that that R is maximum when /3 = — 1. Initially rumor 
size R increases with /3 but after achieving a maximum value for /3 = — 1 it decays 
exponentially. Further, for a = 1, 7 = 2.4 final rumor size R approaches to (beyond 
P = 1.5). Furthermore it is interesting to note that, final rumor size R increased 
with increase of a same as obtained by Eq. ( |T8|) (Fig. |5]). For random inoculation 
g = 0.1, 0.3, 0.5, 0.7, 0.9, the final rumor size R has been plotted against (3 (Fig. |6]). It is 
observed that to get maximum value of R, /3 increases when g increases. Also, maximum 
size of rumor decreases with increase of g. It is because in random inoculation rumor 
threshold value is larger than the threshold in model without inoculation as inferred 
from from Eq. ( |4T1) (Ac > Ac). The sharp decrease in the value of R is seen when rumor 
transmission rate (A) is decrease by 0.5 in comparison to the case where decrease of 
R is shallow when a is decrease by 0.5. Since for A ^ 1 there may be a chance that 
rumor will spread to some extent at any value of g < 1, for very large N. Similar results 
have been observed in the case of targeted inoculation in Fig. [71 But here maximum 
rumor size is much smaller with the inoculation of very less fraction of nodes (e.g. for 
^'=0.25 the final rumor size R is almost suppressed), since rumor threshold in targeted 
inoculation is larger than the random inoculation. For random inoculation strategy, 
the rumor spreading is plotted against time evolution using modified model through 
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-a=0.1, p=-1.1 
- a=03, p=-1 .3 
-a=0.5, p=-1.5 ' 
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(a) a + /3 = -1 



(e) a + /3 = -1 




0.18 
0.16 
0.14 
0.12 
- 0.1 
0.08 
0.06 
0.04 
0.02 





— 1—0=0.1, |3=-0.1 






-e-a=0.3, |3=-0.3 






—^01=0.5, |3=-0.5 






— 01=0.7, p=-0.7 






-B-a=0.9, |3=-0.9 







(b) a + /3 = 



(f ) a + /3 = 




0.18 
0.16 
0.14 
0.12 
- 0.1 
0.08 
0.06 
0.04 
0.02 




-a=0.1, p=:0.9 
-a=0.3, p=0.7 
-a=0.5, p=0.5 
- a=0.7, p=0.3 
-a=0.9, p=0.1 



(c) a + /3 
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(d) a + /3 = 2 



(h) a + /3 = 2 



Figure 2: R(t) vs t (a-d) and S(t) vs t (e-h) for A = 0.8 and different combinations of 
a and /3 parameters 18 



(a) Ac vs a 



(b) Aevs /3 



Figure 3: Threshold (Ac) vs a and /3 for 7 = 2.4 for different size of scale free networks 

simulation results. It is found in Fig. [8]- Fig. [HI for (7 = 0.1, 0.3, 0.5, 0.7 that, if a + /3 
increases from -1 to 2, R will increase since rumor threshold decreases. Further it can 
be observed from Eg. (1321) that, for a + /3 = — 1 and 0, Ac is finite and higher than the 
case where a + /3 = 1 and 2. Therefore, R{t) is almost and grows slowly with time 
when, a + /3 = —1. The growth in R{t) is higher for lower values of but the case is 
reversed for the higher values of g. Similarly in the case of targeted inoculation scheme 
using lower values oi g = 0.05, 0.1, 0.15, 0.2 for a + /3 = — 1 to 2, rumor threshold found 
more than the random inoculation scheme (Eq. (!50|) ) and rumor spreading is supressed 
for inoculation less number of nodes than the random inoculation scheme (Fig. [T2] - Fig. 



8 Conclusion 



In presented study, the modified SIR model has been proposed by considering standard 
SIR rumor spreading model with degree dependent tie strength of nodes and nonlinear 
spread of rumor. The two parameters nonlinear exponent a and degree dependent tie 
strength exponent /3 have been introduced. In modified rumor spreading model, fi- 
nite rumor spreading threshold has been found for finite scale free networks while fixed 
rumor threshold has been found for any size of network when a + /3 + 2 < 7. Ran- 
dom and targeted inoculation schemes have been introduced in the proposed modified 
model. Rumor threshold in targeted inoculation scheme is found to be higher than 
the random inoculation. On the other hand the rumor threshold in random inocula- 
tion is higher than the modified model without inoculation. It has also been observed 
that for scale free networks targeted inoculation scheme is successful in suppressing the 
rumor spreading in the network, since it requires to inoculate less number of nodes 
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(a) 7 = 2.4 



(b) 7 = 3 



Figure 4: Final size of rumor i? vs /3 (A^ = 10^ nodes, a = 0.5 and a = 1) for different 
values of 7 




'0.1 0,2 0.3 0.4 0,5 0,6 0,7 0,8 0,9 1 



Figure 5: Final size of rumor R vs a {N = 10^ nodes, f3 = —1 and A = 1) for 7 = 2.4, 3 

than random inoculation. Further a is found to be more sensitive than /3, as it affects 
more to rumor threshold. Finally it is seen that in real world networks finite rumor 
threshold can be achieved by considering more realistic parameters (degree dependent 
tie strength of nodes and nonlinear spread of rumor). The targeted inoculation scheme 
can be successfully applied to suppress the rumor spreading over scale free networks. 
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Figure 13: R(t) vs t with a + (3 = , \ = 0.6 in targeted inoculation for g= 0.05, 0.1 
(upper) 0.15, 0.2 (lower) 
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Figure 14: R(t) vs t with a + (3 = 1 , \ = 0.6 in targeted inoculation for g= 0.05, 0.1 
(upper) 0.15, 0.2 (lower) 




Figure 15: R(t) vs t with a + (3 = 2 , X = 0.6 in targeted inoculation for g= 0.05, 0.1 
(upper) 0.15, 0.2 (lower) 
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